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We study a one-dimensional model which undergoes a transition between an active and an 
absorbing phase. Monte Carlo simulations supported by some additional arguments prompted as to 
predict the exact location of the critical point and critical exponents in this model. The exponents 
8 = 0.5 and z — 2 follows from random-walk-type arguments. The exponents /3 — u± are found to 
be non-universal and encoded in the singular part of reactivation probability, as recently discussed 
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rf\ • by H. Hinrichsen (|cond- mat/0008179 ) . A related model with quenched randomness is also studied. 
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Q . I. INTRODUCTION 

<D _ 

Recently, non-equilibrium phase transitions have been intensively studied in variety of models. Some attempts 
were also made to classify possible types of phase transitions into universality classes similarly to equilibrium tran- 
3 ■ sitions [00]. However, non-equilibrium transitions are much reacher and puzzling than equilibrium ones and their 
c/3 | understanding is still far from complete. 

A class of nonequlibrium models for which such a categorization seems most feasible are models with absorbing 
states ||. For example, models with unique absorbing states most likely belong to the so-called directed-percolation 
(DP) universality class H. Moreover, models with double (symmetric) absorbing states or with some conservation 
law in their dynamics belong to another universality class, the so-called parity-conserving (PC) universality class [p|. 
Although the above classification is best confirmed for one-dimensional models, it seems to apply to higher-dimensional 
models as well ||. Moreover, models with more than two J7J or even infinitely many absorbing states also fall into 
the above universality classes. (Typically, models with infinitely many absorbing states fall into DP universality class, 
however, certain symmetries ||Q] or conservation laws [[l0| might change the criticality into PC.) 
fvj ' In the present paper we study certain one-dimensional models with infinitely many absorbing states. In the first 

part we study a model defined on a ladder-like lattice (ladder model) and our Monte Carlo simulations show that 
t-H ' this model belongs to neither DP nor PC universality classes. Critical exponents S and z of this model appear to 
be rational numbers (contrary to the DP or PC criticality) and can be obtained using simple random-walk-type 
r ■ arguments. However, the exponents (3 and ^j_(= (3) appears to be non-universal. It turns out that their values are 
encoded in the singular part of the reactivation probability as was recently postulated in the context of related models 
by Hinrichsen flllfl . 

Using certain property of our ladder model we can easily calculate the upper bound on the existence of the absorbing 
phase. Numerical simulations show that this bound most likely give exact location of the critical point in this model. 
In the final part of our paper we examine yet another but closely related model (single-chain model) for which we can 
obtain also the lower bound of the absorbing phase. It turns out that both bounds are the same which exactly locates 
the transition point in this model. The model studied in this part of our paper contains some quenched random 
variables and its critical behaviour is similar to that of the model studied in the main part of our paper. 

o 
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II. LADDER MODEL 

A. Definition and basic properties 

Our model is a certain variant of models recently studied in some other contexts [|9||l2|]. It is defined on a one- 
dimensional ladder-like lattice (see Fig. nk). For each bond between the nearest-neighbouring sites we introduce a bond 
variable w G (—0.5,0.5). Introducing real parameters r and s > 0, we call a given site active when wi'u;2|u>3|' !! < r, 
where w\ and W2 are intra-chain bond variables connected with this site and w^ is the inter-chain bond variable (see 
Fig. |l|a). Otherwise, this site is called non-active. The model is driven by random sequential dynamics and when the 
active site is selected, we assign anew, with uniform probability, three neighbouring bond variables. Nonactive sites 
are not updated. 



First, let us notice that the model has a certain global up-down symmetry. Namely, reversing all bond variables 
w — ► — w we do not change the status of any site, i.e., active (non-active) sites remains active (non-active). The 
same symmetry appears in a closely related model which is defined on a single chain and where activity of a site is 
determined by the product of two intra-chain variables M (in the following we refer to this model as model A). It was 
suggested that this symmetry is responsible for the fact that the model A belongs to the PC universality class H. If 
so, the present model should also belong to this universality class. In the present paper we show, however, that the 
critical behaviour of this model is different. In our opinion, this result does not completely abolish the claims about 
the importance of the up-down symmetry. Namely, as we will show below, the change of the universality class in 
the present model comes from the very intricate mechanism: upon approaching the critical point the system might 
be regarded as composed of two weakly-interacting critical subsystems. We expect that in a generic case such a 
mechanism is absent and models of this kind which are endowed with the up-down symmetry should belong to the 
PC universality class. 

This situation bears some similarity to the equilibrium statistical mechanics, where one expects that short-range 
interacting Ising models will generically belong to the Onsager universality class. It is known, however, that when 
certain additional interactions are included, the critical behaviour of an Ising model, which is then equivalent to the 
eight- vertex model, deviates from the generic case ]13| . 

B. Monte Carlo simulations 

To examine the properties of the above model we used standard Monte Carlo simulations. First, let us describe 
results for s — 2. We calculated the steady-state density of active sites p. The fact that the model is defined on a 
one-dimensional lattice enabled us to examine systems of large size L (up to 5 • 10 5 ). Also the simulation time £mc 
was rather large and typically £mc ~ 10 5 — 10 6 (the unit of time is defined as a single, on average, update per lattice 
site). Our simulations show that for r > the density p remains positive but as soon as r becomes negative the 
system quickly reaches one of the absorbing states and thus p = 0. A simple argument, similar to the one used for 
related models 0], shows that for r < the model gradually generate sites which remains permanently non-active. 
Indeed, let r < and a certain inter-chain bond W3 satisfies the condition 

Kl < ^r/(0.5), (1) 

then for each site connected to W3 the absolute value of the product W1W2W3 is smaller than —r, which implies 
W1W2W3 > r and these sites remain permanently non- active. Since upon updating an active site there is a finite 
probability to satisfy (0), there is a finite rate of creation of such sites (a similar argument can be applied to an 
intra-chain bond). Thus, in agreement with numerical simulations, for r < the system is in the absorbing phase. 
Combined with the computationally observed property that for r > the density p remains positive, it implies that 
r = r c = is the transition point for that model which separates active and absorbing phases. In addition, our 
simulations, which are described below, show that r = is actually a critical point of this model and is accompanied 
by typical power-law characteristics. 

The first evidence of such a characteristic is shown in Fig. 0, which presents the logarithmic plot of p as a function 
of r. The linearity of the plot for approximately three decades confirms the power-law behaviour p ~ r 13 . At the same 
time, it confirms the location of the critical point r c — 0. From the least-square analysis of these data we estimate 
= 0.51(1). 

Another quantity which we measured was the time evolution of the density p for r close to the critical point r = r c . 
One expects that at criticality p ~ t and deviations from the power-law behaviour appear off the critical point. 
The results are presented in Fig. H. For r = the clear power-law behaviour is observed and from the slope of the 
data we estimate S = 0.50(1). 

Other power-law characteristics are obtained from the finite-size analysis. In Fig. Hwe present the size dependence 
of the characteristic time r defined as an average time needed for the system to reach an absorbing state (with a 
random initial configuration). Again, the best linearity is seen for r = and in this case we obtain r ~ L z with 
2 = 2.01(1). Let us notice that this figure limits the allowed values of r c to the very narrow range — 10~ 12 < r c < 10 -8 . 

Finally, Fig g presents the steady-state density p as a function of size L. For r = 10~ 4 and 10~ 6 , p converges to 
the positive value but for r = it decays as p ~ L~P/ UJ - with the exponent /3/v± = 0.99(1). 

For comparison with our model, we quote the values of these exponents for the one-dimensional DP |Q and PC |Fq | 
universality classes: (3 = 0.26486(DP), 0.92(PC); S = 0.159464(DP), 0.286(PC); z = 1.580745(DP), 1.74(PC) and 
z/j_ = 1.096854(DP), 1.83(PC). Thus, our results clearly places the model into a new universality class. 

However, the obtained values of the exponents are not entirely unexpected. In the following we argue that at least 
the values of S, z and of the ratio (3/v± can be inferred from the properties of some other models. First, let us examine 



in more detail the model at the criticality, i.e., at r = 0. Let us notice, that in this case it is only the sign of the 
expression W1W2W3 which determines the state of a site. Since u>§ is always positive, it means that at the criticality 
this term is irrelevant and the state of a site is determined by the product W1W2, i.e., by the product of intra-chain 
variables. The same condition appears in model A and thus the present model at the criticality is equivalent to two 
non-interacting models A. Numerical evidence was already presented M that at r = model A is at the end-point 
of its critical phase which in A model appears in the range < r < r c ~ 0.027. This critical phase is described by 
simple, random-walk related exponents: S — 0.5, z — 2 and /3/v± = 1 Eg ]. On the basis of the above relation and in 
agreement with our simulations we obtain that these random-walk exponents are the critical exponents in our model. 
Let us emphasize that the critical phase in PC models exists in a certain range of a control parameter O] and its 
(random-walk) criticality is described by a different set of exponents than the critical point of such models. In our 
model, there is only an isolated critical point which is described by the random-walk exponents. 

The above relation with model A does not explain the value of the exponent f3, which most likely equals 0.5. This 
is because the relation with A model holds only at r = while the exponent /3 describes the off-critical singularity. 
One can argue, however, that when r is positive but small, then the interaction between these two models A is 
also small. Let us recall now, that for such r model A is critical and thus the present model is equivalent to two 
weakly-interacting critical models A. It is interesting to observe that such an interaction is sufficient to destroy the 
criticality and to keep the system in the active phase with p > 0. Qualitatively, we explain the large (comparing to 
model A) activity as follows: In models with absorbing states a non-active domain can be reactivated only through 
its boundaries. However, the already mentioned weak interaction can reactivate even the interior of such domains, 
which dramatically increases activity of the system. 

C. Reactivation probability 

In a recent paper Hinrichsen has shown that for certain models the singularity of the order parameter is determined 
by singularity of the reactivation probability W ]ll| . His argument relies on the observation that W usually scales 
linearly with the control parameter (i.e., a parameter which measures the distance from a critical point). Thus, in 
models where the density of active sites p scales linearly with W, the singularity of p as a function of the control 
parameter is solely a consequence of the singular behaviour of W. In this subsection we show that Hinrichsen's results 
extends also to the model examined in this paper. 

Let us define the reactivation probability W(r) as a probability that a given active site remains active after an 
update. To compute W(r) we first calculate the probability density P s (z) that wiW2\w 3 \ s — z where w±, W2 and w 3 
are independent and uniformly distributed on (—0.5,0.5). Thus, we have to calculate 



P s (z) = I dwi I dw 2 I -w 3 /s S(z - wiW2W 3 )dw 3 , (2) 



1/2 ,1/2 

dwi I 
-1/2 J-1/2 Jo 

where -w 3 is the probability density of W3 = |w3| s , where w 3 is uniformly distributed on (—1/2,1/2). The 

calculation of the integrals (|2j) is elementary. Pcrfoming integration over w\ , W2 and then over w 3 we obtain 

P s {z) = — / " ^- 2 ln(4|z|M)e(*3~4|z|)^3-^ T {2 s - 1 K4|z|) + .sln(2) + -^] + ^-(4|z|) 1 /-i}, 
s J s — 1 1 — s s — 1 

(3) 

where Q(x) is the unit step function. Having calculated P s (z), the reactivation probability W(r) for positive r is 
given as 

W(r) = f P s (z)dz = 1 + ^-{2 s - 1 [r(ln(4r) - 1) + sr(ln(2) + -?-)] + * ,A ^ S }. (4) 



(l/2)» + 2 



2 



— T {2 s - 1 [r(ln(4r) - 1) + S r(ln(2) + — )] + ^_^(4r) 1 /«}. 



From (|4|) one can see that W(r) is indeed singular at the critical point r — and for s > 1 W(r) ~ r x / s . One can 
see that the numerically found singularity of the order parameter with the exponent /3 = 1/2 for s = 2 corresponds 
to the same singularity of W(r). Thus, the density p in the vicinity of the critical point r — should scale linearly 
with the reactivation probability W(r) and the singular behaviour is only due to the singular behaviour of W(r) as a 
function of r. 

To check the above arguments we performed simulations for s = 4 and 1/2. Let us notice that for s = 1/2 the 
leading term is not r 1 ^ but r ln(4r). In this case (3 should be unity but the logarithmic corrections might substantially 



affect the scaling behaviour. The results are presented in Fig. |[ As estimated from the slope of our data, (3 — 0.26(1) 
for s = 4 and j3 = 1.3(2) for s — 1/2. For s = 4 the exponent (3 is very close to 1/4 which clearly confirms our 
arguments. For s = 0.5 our data are less accurate mainly due to very large relaxation time. Moreover, the logarithmic 
corrections might be responsible for the fact that the 'true' scaling is not clearly seen. 

Let us also notice that at the critical point r — only the sign matters and the behaviour of the ladder model is 
independent on s. Thus, the remaining exponents 5, z and the ratio f3/v± must be the same as in the s — 2 case. 

III. SINGLE-CHAIN MODEL WITH QUENCHED DISORDER 

In this section we examine a simple model which might help us to understand the behaviour of the ladder model 
studied in the previous section. The model is defined on a single chain (see Fig. [j]b). In addition to bond variables 
— 1/2 < w < 1/2 there are on-site quenched variables < v < 1. A given site is defined as active if W1W2V < r, 
where w\ and W2 are variables on bonds connected to this site and v is the corresponding site variable. Once initially 
selected, ^-variables remain unchanged during the evolution of this model (evolving variables are only bond variables) . 
(The already discussed model A corresponds to the case v = 1 on each site.) 

It is elementary to prove the following properties of the this model: 
(i) Similarly to the model studied in the previous section the model for r < generates with a finite rate sites which 
remain permanently non-active. It means that for r < the model is in the absorbing phase. 

(ii) Let r > 0. Let us consider a site for which v < Ar. It is obvious that for any choice of bond variables attached to 
it this site remains permanently active. Since in the thermodynamic limit a finite fraction 4r of sites has v < 4r, we 
obtain that for r > the model is in the active phase. 

As a conclusion, we obtain that r = separates active and absorbing phases of this model. Let us also notice that 
for r = the model is equivalent to the model A (i.e., v variables are irrelevant), which was shown to be critical at 
that point M. Moreover, the fraction 4r of permanently active sites decreases linearly to zero which suggests that in 
this model (3 = 1 g§. 

Monte Carlo simulations of this model strongly suggest that indeed (3 = 1 (see Fig. ||). Similarly to the ladder 
model, at r = only the sign matters and w-variables are irrelevant. Thus, at criticality the behaviour of this model 
is the same as of the A model at r — 0. Numerical simulations of the latter model show Q that in this case z = 2 
and 8 = 1/2 as in the ladder model. 

IV. CONCLUSIONS 

In conclusion, we examined a class of models with infinitely many absorbing states. Critical exponents S, z and 
f3/is± take random-walk values. However, the exponents (3 and v± are non-universal. This non-universality is related 
to the singular behaviour of the reactivation probability W(r). Thus, when as a control parameter of the model we 
choose W(r) rather than r, the order parameter would have the universal exponent (3 = 1. It would be interesting to 
examine this model using, for example, field-theory methods |1!| and to check whether some other universality classes 
among models with absorbing states are possible. 
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FIG. 1. (a) The ladder model. The site (•) is active when u;iu)2|ii>3| s < r. An update of an active site replaces all three 
neighbouring bonds (o) with uniformly distributed random numbers from the interval (-0.5,0.5). (b) The single-chain model. 
The site is active when W1W2V < r and its update replaces only u>i and W2- 

FIG. 2. The log-log plot of the density of active sites p as a function of r. The linear size L in this simulations varied from 
L = 5 • 10 4 to 2 • 10 6 . The solid lines (least-square fit) have slopes corresponding to (3 — 0.26(1) (for s = 4), 0.51(1) (for s = 2), 
1.3(2) (for s = 0.5) and 0.99(2) for the single-chain model. 

FIG. 3. The time evolution of the density p for (from top to bottom) r = 10" 4 , 10" 6 , 0, -10" 10 , -10" 8 , -10" 6 , and 
— 10~ 4 . The simulations were made for s — 2 and L = 10 5 and we checked that for the examined time scale our data are 
basically size-independent. These data limit r c to the range — 10 -10 < r c < 10~ 6 . Even tighter bounds follow from Fig. 0. 

FIG. 4. The size dependence of the characteristic time r for (from top to bottom): r = 10 -6 , 10 -8 , 0, — 10 -12 , — 10~ 8 , 
— 10 -6 , and — 10 -4 . The straight line has a slope corresponding to z = 2.01. Each point is an average of 100 independent runs 
and s = 2. 

FIG. 5. The size dependence of the steady-state density p for r — 10~ 4 (D), 10~ 6 (+) and 0(O). The straight line has a slope 
corresponding to f3/v± = 0.99(2) (s — 2). 
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